Adsorbate-induced substrate relaxation and the adsorbate— adsorbate interaction 
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We formulate the theory of the perturbation caused by an adsorbate upon the substrate lattice in 
terms of a local modification of the interatomic potential energy around the adsorption site, which 
leads to the relaxation of substrate atoms. We apply the approach to CO chemisorption on close- 
packed metal surfaces, and show that the adsorbate-adsorbate interaction and a variety of other 
properties can be well described by a simple model. 



PACS numbers: 68.35.-p, 68.35.Md, 82.65.My 

Several direct and indirect (through-the-substrate) 
mechanisms can lead to an effective interaction between 
adsorbates on metal surfaces. The possible role of 
adsorbate-induced substrate relaxation was considered 
some time ago by Lau and Kohn using an elastic 

continuum model of the surface. They concluded that 
the resulting interaction was repulsive between identi- 
cal adsorbates, varied as p^^ with separation p, and was 
inversely proportional to the shear modulus of the sub- 
strate. For adsorbates separated by a few atomic spac- 
ings, however, the continuum theory may not be valid. 
More recently, in a series of papers Kevan et al. have 
determined adsorbate-adsorbate interaction energies for 
CO at several metal surfaces, using a transfer-matrix 
analysis of thermal desorption spectra. They have also 
qualitatively discussed the adsorbate-induced strain as a 
possible mechanism of the adsorbate-adsorbate interac- 
tion at intermediate range, i.e., several substrate atoms 
apart. 

The potential energy of the atomic lattice of a solid in 
the harmonic approximation can be written as 
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where x^{i) is the fi-th component of the displacement 
of the z-th atom from the equilibrium position. Now as- 
sume that an impurity is created by replacing one atom 
with a different species. (In the general discussion we 
talk about an "impurity" , although we are primarily in- 
terested in the chemisorption case, where the adsorbate 
also introduces additional degrees of freedom. The gen- 
eralization to the latter case is straightforward.) The 
potential energy after the replacement can again be ex- 
pressed in the form (^ , but the new equlibrium positions 
x'^(i) are in general different: 
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The new dynamical matrix D'^j^{i,j) is also different, and 
a constant term appears which shifts the energy mini- 
mum. Using (|^), the potential energy of the system after 
the impurity has been introduced can be expressed in the 
coordinates x^ii): 
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i.e., the dynamical matrix is modified, and linear (force) 
and constant (energy shift) terms appear. 

The foregoing considerations depend only upon the as- 
sumed stability of the solid, i.e., the existence of the 
minimum of the potential energy. For our application 
to chemisorption, we specifically assume the following 
properties: (a) the force terms ^)i(j) are nonzero only 
for a small number of substrate atoms around the ad- 
sorption site; (b) similarly, only a few elements of the 
dynamical matrix Z?p,y(i,j) change, if any; (c) the effect 
is linear, so that the chemisorption of another molecule 
(and, consequently, of a third, a fourth, etc.) can be de- 
scribed by the same set of parameters, of course centered 
around the new adsorption site. The condition (c) is 
rather restrictive and excludes systems which reconstruct 
at large adsorbate coverage, as well as those where other 
interactions are important, such as the direct adsorbate- 
adsorbate repulsion, the electrostatic dipole-dipole in- 
teraction, or the "chemical" competition for the same 
electronic orbitals in the substrate. Estimates, however, 
show that such interactions are usually weak and short- 
ranged, while we are interested in medium-range interac- 
tions (second nearest neighbor and beyond) of nonionic 
adsorbates. 

In this paper we take into account the in-plane re- 
laxation within the first atomic layer of the substrate 
induced by a chemisorbed species. The relaxation in 
the perpendicular direction can, of course, be equally 
strong, but it does not contribute much to the effective 
adsorbate-adsorbate interaction and to the adsorbate- 
induced surface stress. Also, we do not consider the cou- 
pling to internal adsorbate coordinates, which was dis- 
cussed in Ref. 1^ in an application to the damping of 
adsorbate vibrations. 

For a hexagonal layer of atoms, we write the potential 
energy as 
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where rj = (xi,yi) is the in-plane displacement from the 
equlibrium position of the i-th atom. The term Ki de- 
scribes a central atom-atom interaction, and the term 
K2 binds atoms to their equilibrium positions, simulat- 
ing the interaction to lower atomic layers. Without it, the 
model would be too "soft" to long-wavelength perturba- 
tions. The trade-off is that the lowest phonon frequency 
becomes finite, i.e., there are no true "acoustic" modes, 
but this has little influence on the relaxation energy and 
other static quantities calculated in this work. 

Now assume that an atom or a molecule chemisorbs on 
top of the atom i = 0, and that the induced change of 
the potential energy is 



r = ^iAXi[fo,.(ro-r,)p 
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where the first term describes the change of the force con- 
stant between the atom and the six surrounding atoms, 
and the second is a linear force term. In the following we 
drop the AKi and Vq terms, which are beyond the level 
of accuracy assumed in this work, although, in general, 
the AKi term ought to be included. 

The definition of the chemisorption-induced surface 
stress T is PI 
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where A is the surface area, Scx the strain, and SW the 
difference of the work involved in straining a clean surface 
and a surface with adsorbates. In our model, we obtain 
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where 9 is the adsorbate coverage. (If AKi is different 
from zero, an additional factor of order unity appears in 
the above expression.) 

Irrespective of the sign of the force term k2/a, there 
is always an energy gain due to the relaxation of the 
surrounding atoms, i.e., the minimum oi V + V is less 
than zero. If another molecule adsorbs on a nearby site, 
the forces induced by the two adsorbates act in oppo- 
site directions and the relaxation is less complete than 
with adsorbates far apart, which leads to an effective in- 
teraction. We calculate the interaction energies by com- 
paring the minimum of the potential energy V + V' for 
a single adsorbed molecule with the minimum for two 
molecules adsorbed on, in turn, second, third, fourth, 
and fifth nearest-neighbor sites, see Fig. ||(a). (Through- 
out the paper we assume that there is a strong repul- 
sion between first nearest neighbor adsorbates caused by 
"chemical" effects, and do not consider them.) 




FIG. 1. (a) Chemisorption of CO, black circle, into an 
on-top site on Pt(lll). The chemisorbate causes a relax- 
ation of the adjacent Pt atoms, denoted by arrows, which 
propagates to further first-layer atoms via elastic forces. The 
neighboring adsorption sites are indicated by numbers, (b) 
The regular (^/3 x \/3)R30° structure at a coverage 6 = 0.33. 
Owing to symmetry, no relaxation of substrate atoms along 
the surface plane is possible. 



We first consider the chemisorption of CO on a Pt(lll) 
surface. There is a large amount of tensile stress in the 
first atomic layer of the clean Pt(lll) surface pO| , pl| . 
Although the surface does not reconstruct at room tem- 
perature, a reconstruction is observed at high temper- 
atures in the presence of saturated Pt vapor ||T^. CO 
adsorbs initially into on-top sites of Pt(III), but the 
energy difference for the adsorption into bridge sites is 
obviously small, since some bridge adsorbates are found 
already at coverages 9 above 0.15 jl^. Some authors 
report (Ref. [|l^ and references therein) that a regular 
(V3xy3)R30° structure at a coverage 6 = 1/3 is formed. 
Fig. 0(b), but others claim that the densest structure of 
exclusively on-top adsorbates exists at 6* = 0.29 and 
that further chemisorption occurs into bridge sites. The 
regular structure aX 9 — 0.5 contains an equal number of 
on-top and bridge adsorbates. 

We have chosen the values of the parameters Ki, K2, 
and ^2/0; which give good agreement with experimen- 
tal data on adsorbate-induced surface stress |^ and with 
low-coverage interaction energies , as shown in Table |. 
We had to choose a rather small value for the force con- 
stant Ki between atoms in the first surface layer. (The 
value of K2 has little effect on the results. The much 
larger value of Ki used in a similar model in Ref. 
was an overestimate.) The reduction from bulk values 
is characteristic of many close-packed noble-metal sur- 
faces but the reduction we find is larger than sug- 
gested in the surface phonon calculation in Ref. [ p^ . 
Consequently, the highest resulting vibrational frequency 
of two-dimensional phonons of the first surface layer is 
about 5 meV, by about a factor of two smaller than the 
frequency of surface phonons along the edge of the Bril- 
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louin zone of around 10 meV calculated in Ref. As 
discussed above, it is possible that part of the softening is 
localized around the adsorption site only, the term AKi 
in Eq. (|^). The values in Table || show that the repul- 

TABLE I. Interaction energies W„nn (in K) and the in- 
duced surface stress r (in N/m) for on-top adsorption of 
CO on Pt(lll). W„NN = 2Eo — -EnNN, where Eq is the 
relaxation energy of a single CO molecule, EnNN is the re- 
laxation energy for two CO molecules adsorbed at the n-th 
nearest-neighbor sites, etc. Theoretical values are calculated 
using Ki=4 N/m, K2 = 2.5 N/m, fe/a = 0.3 x 10"^ N. The 
nearest-neighbor distance is a = 2.76 A. 
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''Ref. M. low CO coverage. 
''Ref. Ij, CO coverage of 0.33. 
'^Ref. PI, CO coverage of 0.33. 



sive interaction is strong between CO adsorbed on sites 
lying along rows of substrate atoms, and weaker for ad- 
sorbates separated by hollows, even if they are less far 
apart. In our opinion, the rather large interaction energy 
between fourth nearest neighbor adsorbates in Ref. [Q is 
influenced by the contributions from more distant sites, 
which were not included in their analysis. 

Clean nickel (111) surfaces do not reconstruct. Un- 
like some earlier claims, it is now accepted that at low 
temperature CO chemisorbs initially into threefold hol- 
low sites 0,|8|. At room temperature, some bridge 
and on-top sites seem to be occupied even at low cov- 
erages M. At a coverage 6 — 0.33, CO forms a reg- 
ular (V3x \/3)R30° structure jl^, but it is not clear 
whether the molecules adsorb into fee or hep positions. 
At 6 — 0.5, a regular c(4x2)-2CO structure is formed, in 
which an equal number of fee and hep sites is occupied 
[|9|-|l|. Fig. |(b). The top layer of Ni atoms shows buck- 
ling in that nonequivalent atoms have different vertical 
relaxation. The CO molecules are slightly tilted from the 
direction perpendicular to the surface. The initial heat of 
adsorption at room temparature is 130 kJ/mol It de- 
creases slowly at first, to around 122 kJ/mol dX — 0.33 
and 112 kJ/mol at 6* = 0.5. However, the fact that not all 
adsorption sites are equivalent, and the rather large stan- 
dard deviation of experimental data make the interpreta- 
tion of the coverage dependence of the heat of adsorption 
uncertain. 

We describe the first layer of Ni atoms by the same 
potential as for Pt, Eq. (|^). The interaction terms are 
similar to Eq. (||), but the adsorbate is in a threefold 
hollow site and the sums run over the three surrounding 
Ni atoms (Fig. H). The induced surface stress is 




FIG. 2. (a) Chemisorption of CO, black circle, into a 
three-fold hollow site on Ni(lll). The relaxation of the ad- 
jacent Ni atoms is denoted by arrows. The numbers indicate 
possible adsorption sites of the same kind, either fee or hep. 
(b) The regular c(4 x 4)-2CO structure on Ni(lll) at a cov- 
erage 9 = 0.5. An elementary cell (grey area) contains one 
adsorbate in the fee and one in the hep position, d and e. 
The relaxation of Ni atoms in intermediate rows, a and c, is 
denoted by arrows. 

where the nearest-neighbor Ni-Ni distance is a = 2.49 A. 
From the experimental data r = —0.55 N/m at 6* = 
0.33 ^ we have estimated k2/a = 1 x 10"^'' N. In a 
LEED analysis of the c(4x 2)-2CO structure which forms 
at 9 = 0.5, Mapledoram et al. found that the lateral 
displacement of the Ni atoms next to an adsorbate, a and 
c in Fig. |^(b), was 0.03 A. We have reproduced this value 
by taking Ki = 6 N/m and K2 = 2 N/m. (As before, the 
value of K2 is of lesser importance.) This is a consider- 
able reduction from the bulk values, but not as large as 
for the Pt(lll) surface, in agreement with the fact that 
the Ni(lll) surface seems less prone to reconstruct than 
Pt(lll). The energy gain per adsorbate is 22 K. 

Using these values of the parameters, the relaxation 
energy for a single adsorbate is only 42 K, which is 30 
times smaller than that we have found for CO/Pt(lll). 
(This reduction can be well estimated treating each sur- 
rounding Ni atom as an independent oscillator: The force 
^2/0: is three times smaller than in the CO/Pt(lll) case, 
the Ni-Ni force constant is 50% larger, and there are 3 
surrounding atoms instead of 6.) Furthermore, the dis- 
placements of Ni atoms around the hollow adsorption 
site are not along chains of atoms, and do not cause a 
large displacement of other atoms. The calculated effec- 
tive interaction energies for second nearest neighbors and 
beyond are therefore only a few K, too small to be ob- 
servable. In our opinion, the interaction energy of 100 K 
between second neighbor adsorbates suggested by Skel- 
ton et al. ||] is either a combined result of other mecha- 
nisms or an artefact of their procedure. We note that an 
earlier study reported that there was essentially no 
interaction already between second neighbor adsorbates. 
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Interaction energies have also been determined for CO 
chemisorbed on Rh [0|2^ and Cu ^ surfaces. We discuss 
these systems only qualitatively, since the proposed val- 
ues are still uncertain, and there is no quantitative data 
on other adsorbate-induced properties. Wei et al. ||^ es- 
timated that W2 = -100 K and W3 = 150 K for the 
on-top chemisorbed CO on Rh(lll). An earlier mea- 
surement by Payne et al. |25| reported W2 = 170 K and 
W3 — —85 K. In our model, the relative magnitudes of in- 
teraction energies for on-top adsorbates on fee (111) sur- 
faces are always similar to those found for CO/Pt(lll). 
In particular, we expect large repulsion between third 
nearest neighbor adsorbates which lie along a chain of 
substrate atoms. In this respect the values proposed in 
Ref. seem more probable, although the origin of the 
attractive W2 (if it is real) is not clear. For the on-top 
CO on Cu(lll), the same authors found W2 = 107 K, 
> 800 K, W4 = 155 K §. The value of W3 seems 
too large compared with the other two, but otherwise the 
results are quite similar to CO/Pt(lll). 

The square lattice of the first layer of the (100) surface 
of the fee bulk is not well described by a purely pairwise 
potential similar to Eq. (^), since in the absence of angu- 
lar force constants and of coupling to lower layers, only 
the ad hoc term K2 ensures the stability. Nevertheless, 
the trends in the interaction energies between on-top ad- 
sorbates can be deduced by analogy with (111) surfaces. 
We expect a repulsive interaction between third near- 
est neighbor adsorbates which lie on the same chain of 
atoms, and no interaction between second nearest neigh- 
bors which lie diagonally on different chains, because the 
forces on adjacent substrate atoms are orthogonal. In the 
latter case, even a weak attraction is possible owing to 
partly coUinear displacements induced by the two adsor- 
bates on more distant surface atoms. Indeed, the values 
W2 = and W3 = 400 K were found for Rh(lOO) §, 
while the values suggested for Cu(lOO) were W2 = —33 K 
and 14^3 = 13 K |6|. It is interesting that the contin- 
uum elastic theory |^] also gives a strong repulsion in 
the (110) direction and possibly a weak attraction in the 
(100) direction between adsorbates on (100) surfaces of 
noble metals. 

We have shown that the adsorbate-induced substrate 
relaxation leads to an adsorbate-adsorbate interaction 
which is quite long-ranged and has a nonmonotonic 
distance dependence, being particularly large between 
molecules adsorbed in on-top positions along a chain of 
surface atoms. The same mechanism also leads to other 
observable effects, such as the adsorbate-induced sur- 
face stress and the relaxation displacement of substrate 
atoms. Quantitative agreement with experiment can be 
obtained for CO adsorbed on Pt(lll) and some other sur- 
faces using a simple model, assuming that force constants 
between atoms in the first surface layer are considerably 
weaker than in the bulk, which is a known property of 
many close-packed noble metal surfaces. The results em- 



phasize the importance of allowing the full substrate re- 
laxation in the first-principle calculations of chemisorp- 
tion on metal surfaces. 
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